(N 
O 
O 
(N 

> 

o 
O 

m 



Experimental study of the spatial distribution of quantum correlations 
in a confocal Optical Parametric Oscillator 

M. Martinelli, N. Treps, S. Ducci, S. Gigan, A. Maitre, C. Fabre 
Labomtoire Kastler Brossel, Universite Pierre et Marie Curie, 
case 74, 75252 Paris cedex 05, France. 
(Dated: February 1, 2008) 

Abstract 

We study experimentally the spatial distribution of quantum noise in the twin beams produced by a type 
II Optical Parametric Oscillator operating in a confocal cavity above threshold. The measured intensity 
correlations are at the same time below the standard quantum limit and not uniformly distributed inside the 
beams. We show that this feature is an unambiguous evidence for the multimode and nonclassical character 
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I. INTRODUCTION 



To date, almost all the experiments investigating the quantum properties of the states of light 
produced by Optical Parametric Oscillators (OPO) have been performed through the measurement 
of the total intensity of the generated fields, obtained by integrating on the detector the whole 
wavefront intensity distribution. Such measurements have put in evidence the squeezed vacuum 
character of the output of a degenerate OPO below threshold |1|], the quantum intensity correlations 
between the signal and idler beams (twin beams) 0], and the bright squeezing of the output of an 
optical parametric amplifier (OPA) |^ and of the pump beam reflected by an OPO [Q. 

In the last years, the interest has turned to the spatial aspects of quantum fluctuations, in par- 
ticular because they open new perspectives in the quantum information field : spatial features offer 
the possibility of parallel processing and multichannel operation, each part in a transverse section of 
a beam playing the role of a channel. The concepts of temporal squeezing and correlations of light 
beams as a whole have been extended to spatio-temporal properties for the local quantum fluctua- 
tions in the plane perpendicular to the propagation : squeezing of the temporal fluctuations inside 
a small part of the transverse plane, or temporal quantum correlations between different transverse 
areas of the beam are some of the manifestations of the transverse quantum properties of light. It 
has been shown that, in order to obtain strong spatial quantum effects, one needs to use multimode 
nonclassical states of light, in contrast to the single mode operation of the experiments quoted 
above. It has also been shown that such nonclassical multimode states can be used to improve the 
optical resolution and to measure small transverse displacements of a light beam below the stan- 
dard quantum limit [^, 0. 

Our investigations of the spatial distribution of multimode states of light produced by OPOs are 
related to the pioneer studies of the stochastic spatial noise performed in gas lasers P], and more 
recently in semiconductor lasers (diode lasers and VCSELs) 0, for which the spatial distribution 
of the quantum fluctuations of the output beam has been studied and interpreted in terms of a 
superposition of several Hermite-Gauss modes. They must also be connected to the studies of spa- 
tial quantum effects for the photons generated by spontaneous parametric down-conversion, at the 



photon counting level, for example to the recent demonstration of spatial antibunching [p!0 



Parametric interaction in a non-linear crystal has been theoretically shown to be a very efficient 
source of multimode nonclassical states of light. Many theoretical studies |Tl|, have shown in 
particular that Optical Parametric Amplifiers (OPA) generate multimode squeezed states. OPOs 
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have also been shown to be a source of multimode squeezed states provided that they use optical 
cavities with degenerate transverse modes : in particular, detailed theoretical investigations have 



been performed on sub-threshold OPOs in planar |T3, WM, or quasi-planar cavities p!5|, [16 



From an experimental point of view, such cavities, which are close to instability, are quite difficult to 
handle and have very high oscillation thresholds. Confocal cavities are much easier to operate, and 
still exhibit interesting nonclassical transverse effects, such as multimode squeezing in the degenerate 
case below threshold The multimode transverse operation of confocal OPOs above threshold has 



been theoretically investigated at the classical level [T^, and experimentally put in evidence in [|T^], 
where it has been shown that the field emitted by a confocal OPO above threshold can be described 
as a superposition of a great number of TEMpg modes. The quantum description of this regime 
is unfortunately a difficult task, as the device switches from an uncoupled regime of the different 
transverse modes to a strongly coupled regime when passing from the sub-threshold emission to the 
regime of intense output. 

This paper is devoted to the study of the spatial distribution of quantum fluctuations of such a device, 
and shows that it generates multimode nonclassical states of light. In section II, we give the precise 
definition of multimode quantum states of light. We present in section III a criterion unambiguously 
characterizing such states in our experimental configuration, namely the measurement of the spatial 
distribution of the intensity correlations between the signal and idler beams. Section IV presents the 
experimental setup and the procedure for measuring these correlations. Finally, section V presents 
the experimental results proving that the signal and idler beams emitted by the confocal OPO are 
multimode and spatially quantum correlated. 



II. INTRINSIC CHARACTERIZATION OF MULTIMODE AND SINGLE MODE BEAMS 
A. Definition of single mode and multimode beams 

In the studies on optical patterns, a beam is said to be multimode (in the TEMpg basis), when 
the far field patterns appear to be different from the near field patterns. Actually, this feature is a 
proof that the beam consists of a superposition of several TEMpg modes. The electric field positive 
frequency envelope of such a light beam, E^^\f, z), expanded in the basis of TEMpg modes Up^g{f, z), 
writes : 

E^+\f,z) = ^ap,gUp,g{f,z). (1) 

P.9 
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with more than one ap^q coefficient different from zero. Nevertheless, if the ap^q coefficients are fixed 
(i.e. if we deal with a coherent superposition of modes and not a statistical one), one can always 
define a new mode Vo{r, z) : 



and construct a new orthonormal basis of modes {fi(r, z)} in which is the ffist element. In this 
new basis, the field which appeared multimode in the Gauss-Laguerre basis, is proportional to vq 
and is then single mode. This simple reasoning at the classical level seems to show that the single 
or multimode character of a beam having a well defined and fixed amplitude distribution depends 
on the choice of the basis, and is not an intrinsic property. A TEMqq laser beam is single mode in 
the Laguerre- Gauss basis, and multimode in the basis of transverse plane waves. We will show here 
that this statement is no longer true when one describes the beam at the quantum level. 
We will now define a single mode quantum state of light in the following way : | > is a single 
mode quantum state of light if there exists a basis of modes {fi(r, z)} in which it can be written as 



where the first transverse mode vq, whatever its shape, is a non-vacuum state | \E'o >? and all the 
other modes are in the vacuum state. We will call intrinsic multimode states all states defined 
by a ket | ^ > which cannot be written as @ in any basis. We will give in the next subsection a 
characterization of such single mode or multimode states which is independent of the basis used to 
describe it. 

This striking difference between the classical and quantum description of a multimode state comes 
from the fact that a quantum state gives information on the mean electric field, but also on the 
statistical distribution of quantum fluctuations : a single mode beam is characterized by a well 
defined transverse variation which carries all its transverse information. Consequently, in such a state 
the transverse distribution of quantum fluctuations can be deduced from the transverse variation of 
the mode itself. In contrast, in intrinsic multimode states, the spatial distribution of fluctuations 
and correlations cannot be deduced from the structure of the mean field. 

Let us finally mention that this problem is modified if one considers stochastic fields, i.e. having 
classical fluctuations. In this case also, the comparison between the spatial (or temporal) distribution 
of the mean field and of the classical fluctuations is still an intrinsic tool to determine whether the 
field can be considered as single mode or multimode |0. We will not consider further this problem 




(2) 



^ >=| ^o>®|0>(8)|0> (g)... 



(3) 



here. 
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B. Characterization of a single mode quantum state 



Let us consider a single mode state of light | \1/ >, written in the adapted basis {vi} as (H), and 
let us call di the annihilation operator of photons in the mode Wj. Let us now introduce any other 
mode basis of the transverse plane {wj}, and the corresponding annihilation operators bj. There 
is a unitary transformation relating the two basis, and correspondingly the two sets of annihilation 
operators : 

j 

From equations (|^) and one deduces that : 

>= U,odo I ^'o > (5) 

We thus obtain a specific property of a single mode quantum state : the action on it of all the 
annihilation operators of any given basis of transverse modes gives vectors which are all proportional. 
One can show that this feature is a necessary and sufficient condition to be a single mode state . 
In contrast, for an intrinsic multimode state, there exists at least one basis of modes in which this 
property is not true. The condition for single mode states being very restrictive, the set of single 
mode states is a very small subset of the general Hilbert space. The complementary of this subspace, 
namely the set of intrinsic multimode states, is therefore much larger. 

This characterization of single and multimode states appears to be quite mathematical, and seems 
difficult to implement experimentally. We give in the following a more convenient property of single 
mode states, which is not a necessary and sufficient condition, but which can be submitted to an 
experimental check with our set-up. 



III. CHARACTERIZATION OF SINGLE MODE AND MULTIMODE TWIN BEAMS 

A. Partial measurement of intensity fluctuations on a single mode beam 

Let us consider an intensity, or photon number, measurement using a detector of variable trans- 
verse area Sa- We have shown in reference Q that if the light is in a single mode state, the variance 
of the photon number fluctuations measured with this partial detector is given by : 

ANl =< Na > < N,^, >) (6) 
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where < Na > is the mean number of photons detected on the area Sa, < Ntot > the mean number 
of photons detected on the whole transverse plane, and ANtot the corresponding variance. One sees 
that the intensity noise normalized to shot noise, varies linearly with the quantity T = -^^f^- 

This formula can be understood by considering that in a single mode state the photons are randomly 
distributed in the transverse plane, so that a partial detection will introduce sorting noise in the 
detection of the photon number, exactly like when one introduces a linear loss of value T in front of 
the detector. This intuitive picture is no longer true for multimode beams. 



B. Partial measurement of intensity correlations on single mode twin beams 

In a non degenerate OPO, the emitted signal and idler fields present strong intensity correlations at 
the quantum level ("twin beams"). In particular, the variance of the difference Ndif^tot = ^i,tot — ^2,tot 
between the signal (labeled 1) and idler (labeled 2) intensities detected over the whole transverse 
plane is smaller than < Ni^tot > + < N2^tot >, which is the shot noise for the sum of the signal and 
idler beams. When the beams produced by the OPO are both single mode, one can easily show from 
1^ that the following formula, similar to (|^), holds in case of a partial photodetection of signal and 
idler beams by two detectors having the same areas and positions Sa in the two beams : 

^Nlf^A = (< N,,A > + < N,,A >)+7^#^^^^^#^^(A^i/,to -(< N,,ot > + < iV2,to* >)) 

(7) 

Similarly to the interpretation of formula (^), this formula can be understood by considering that in 
single mode twin beams there is no spatial correlations between the twin photons inside the signal 
and idler beams, the photons being randomly distributed inside the two beams. If, instead of small 
detectors, we use detectors with a broad enough area to detect the whole beams, but preceded 
by a diaphragm, or iris, of variable transmission T, expression (|^ shows that in single mode twin 
beams the relative noise on the intensity difference, ua = <jVi a>+<'a^2 a> ' ^ linear function of the 
transmitted mean intensity. In contrast, a nonlinear variation of this quantity will be a signature for 
intrinsic multimode twin beams, for which the twin photons are not randomly distributed inside the 
two beams. 
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C. Partial measurement of intensity correlations on single mode twin beams undergoing 
different losses 



We have assumed so far that the detected mean intensities of the output beams are always equal, 
even in partial photodetection. In real experiments, this is not exactly the case, and imbalances of 
several percents are commonly measured, especially in confocal OPOs, where the patterns observed 
in the signal and idler beams are different. The imbalance will affect the intensity fluctuations of the 
two beams in a different way. The imbalance of the whole beams is due to the fact that the signal 
and idler beams, produced in exact equal amounts in the nonlinear crystal, experience different losses 
in their propagation. Furthermore, in the case of a partial photodetection, one must also take into 
account that the signal and idler beams may have different shapes. It is important to precise the 
way to define the "twin character" of such imbalanced beams in a partial measurement. We will first 
introduce the normalized intensity difference noise n when the whole beams are measured (Fig.[l|a). 
We will then consider a partial measurement (Fig.|I|b) and define the direct normalized intensity 
difference noise n^, with no corrections, as it can be directly calculated from the measurement. We 
will finally define the corrected normalized intensity difference noise Ucorri an attempt to recover the 
value of n from assuming that the beams are single mode. 

Assuming that the mean signal and idler fields are real, we can write the field intensity operator, 
in the small fluctuation approximation, as : 

5im = 2Em5E^ me {1,2} (8) 

where SEm is the quadrature fluctuation operator in the amplitude direction. We now introduce the 
normalized quantity n, which is the variance of the intensity difference fluctuations normalized to 
the total shot noise, and can be written as : 

4(^1 + ^,)^^E^ 

where A^i?^ is the noise variance of the vacuum field, n > 1 for beams with classical correlations 
(obtained for example by using beamsplitters) and n < 1 for quantum correlated beams ("twin 
beams"). Let us now assume that the two beams are subject to losses before being incident on the 
photodetectors (Fig.pD). We call E°^^ the fields which are incident on the photodetectors and i^* 
their intensities, and the normalized variance of their intensity difference : 




5E,2{t) 

FIG. 1: measurement of twin beams correlation in the imbalanced case, z^"*: direct difference measurement, 
difference measurement after gain corrections. 

In the transverse single mode case, losses of any origin and partial photodetection have the same effect 
on the fluctuations mechanism, identical to the effect of a beam splitter of transmission (m=l,2) 
(with = i'^^/im)- The field intensity measured by the photodetectors present fluctuations 
5i^*(t) given by: 

5C\t) = 2tlEm5Em + 2t™V(l - tDEmSKrr. mG {1,2} (11) 

Knowing the effective loss coefficient t^, it is possible to adjust the gain on the acquisition 
channels at the value 1/t^, as shown in fig. ^ so that the mean values of the corrected intensities 
9m{i°m) equal to their value without losses {im)- If one assumes that the two beams are single 
transverse mode, it is possible to recover the original normalized intensity difference noise n of eq. 
(P) by subtracting the effect of the vacuum fluctuations SEym- Writing A^z°"* the fluctuations of the 
difference between the two acquisition channels : 

AH";' = {{g,5lT' - (12) 
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One easily shows that the quantity 

is equal to n in the single transverse mode case. The two last terms of (|T3|) correspond to the 
corrections arising from vacuum fluctuations introduced by the losses or the partial detection. 



D. Multimode twin beams 



Since for single mode signal and idler b constant value ricorr with the transmit- 

tance tm for a partial measurement (or a linear variation of n^) is a good indication (but not a proof) 
of the single mode character of the signal and idler mode beams. In contrast a non constant value of 
Ucorr-, or a nonlinear variation of n^, with respect to the losses is a an unambiguous signature of the 
multimode character of the two beams generated by the OPO. Furthermore, if these quantities are 
smaller than 1, we can conclude that we are in presence of multimode nonclassical beams. Strictly 
speaking, as the vacuum correction in ricorr has been derived assuming single mode beams, ricorr do 
not anymore correspond exactly to a noise correlation in the case of multimode beams. Moreover, 
for a partial measurement, the mean intensity of both beams may be different, and neither rid nor 
ricorr are perfectly suited for the exact characterization of quantum intensity correlation. For these 
reasons, in the case of a partial measurement of a multimode beam, more sophisticated criteria 
should be investigated, that will be developed in a forthcoming publication. 



IV. EXPERIMENTAL SETUP 



For OPOs operating in cavities with degenerate transverse modes, the formation of complex 
spatial structures has been theoretically [^1], and experimentally [0, ^ studied. Furthermore, 
quantum multimode operation of OPO has been already theoretically predicted [Q. But, up to 
now, to our knowledge, there has been no experimental demonstration that the emission of an OPO 
is intrinsically multimode. The experiment presented here investigates the spatial distribution of the 
intensity correlation between the signal and idler beams emitted by a confocal OPO above threshold 
in order to identify their single mode or multimode character. The setup, shown in Fig.^ can be 
divided in two parts: the triply resonant OPO and the acquisition system of the intensity fluctuations 
of the beams. 
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FIG. 2: Experimental setup 

A. The triply resonant OPO 

In order to match the desired requirements of power and stability of the experiment, a dedicated 
setup was made for the generation of the pump beam for the OPO. We used an ultrastable single 
mode Nd:YAG (yttrium aluminum garnet) laser with 350 mW of output power to lock a flash-lamp- 
pumped Nd:YAG laser operating in a ring cavity. This configuration p9| , p3| gives a stable single 
mode beam with an output power of 3.5 W at 1064 nm. This beam is injected in a semi-monolithic 
frequency doubler, using a MgOiLiNbOs crystal and a concave mirror. The output of the second 
harmonic generator reproduces the frequency stability of the injected Nd:YAG laser and gives a total 
output power of 1.3 W at 532 nm. The pump is then injected inside the OPO with a waist equal to 
twice the TEMqq waist of the OPO cavity mode like in reference |T^. 

The nonlinear crystal used inside the OPO is a walk-off compensated KTP (potassium titanyl phos- 
phate) cut for type II phase matching. It is made by two 5 mm long crystals that are optically 
contacted, with their orientations compensating the walk-off effect. The triply resonant OPO uses 
two plane-concave mirrors, whose curvature is R = 100 mm. The input coupler (IC) has a transmit- 
tance of 10% for the pump (532 nm) and high reflectance for the Nd:YAG wavelength (1064 nm). The 
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output coupler (OC) has a high reflectance at 532 nm and 1% transmittance at 1064 nm. Measured 
finesse value of the cavity with the crystal is 40 at 532 nm and 300 at 1064 nm and temperature is 
controlled such that signal (o-polarization) and idler (e-polarization) emitted by the OPO are close 
to frequency degeneracy. The cavity length can be tuned over a few free spectral ranges with a 
piezo-electric ceramic (PZT) attached to the input coupler. A coarse control of the length of the 
OPO is performed by means of translation stages fixed on both cavity mirrors. A fraction of the 
infrared output of the Nd:YAG laser can also be injected into the cavity for alignment purposes and 
to check the transverse mode separation. The OPO works close to the confocal configuration. In 
order to define the region of confocality, it is useful to remember that the distance L between two 
spherical mirrors for a confocal cavity depends not only on the ray of the concave surface but also 
on the refractive index of the medium inside the cavity. Therefore, considering the diffraction effects 
inside a crystal of length £ and refractive index n in a cavity made by two spherical mirrors of ray 
R, the distance L between the mirrors for the confocality will be 



L,onf = R + (il-l/n) (14) 

Since the refractive index for pump, signal and idler are different (respectively 1.7881, 1.8296 and 
1.7467 [^), the length of the cavity L^onf for which the exact confocality is obtained, is different for 
each of the three modes (104.41, 104.28 and 104.53 mm). Transverse degeneracy is achieved when 
different transverse modes are resonant for the same cavity length. Due to the width of the resonance 
peaks, transverse degeneracy can be obtained even if the cavity is not exactly confocal. The range 
of confocality can be defined as the region of the cavity length where the fundamental and the first 



transverse even mode separation is smaller than the cavity bandwidth, and, following [^, can be 
expressed as 

ttR 

|L-Lco„/|<— (15) 

where F is the finesse of the cavity. In our case, the range of confocality for signal, idler (and pump) 
is larger than the difference between the confocal lengths for the signal, idler and pump modes. 
Transverse degeneracy, essential for multimode emission, is then achievable at the same time for the 
three modes. Taking the average of the confocal length for signal and idler, we will consider the 
confocal length as Lconf = 104.41mm, and express the cavity length in terms of AL = L — Lconf- 
Although the threshold for oscillation is quite low and around 30mW, significant thermal effects take 
place inside the cavity even close to threshold, because of the non-negligible absorption of the green 
light(3% at 532nm) and of infrared light (0,45% at 1064nm) : thermal lensing changes the total 
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Gouy phase shift added to the wave in a round trip p6[ and even a sub-confocal cavity can become 
transverse degenerate if the pump power injected into it is sufficiently high. Therefore the confocahty 
iQf is obtained for cavity lengths which are shorter than the confocal length Lconf, defined above for 



a "cold" OPO with no thermal effects. 



B. Detection and acquisition 

The details of the setup used for studying the transverse distribution of noise in twin beams 
cavity can be seen in Fig.|^. Near the output coupler, a dichroic mirror (DM1) is used to eliminate 
the residual green light transmitted through the OPO cavity. A small fraction (2%) of the infrared 
light that is reflected by this mirror and transmitted by a second dichroic mirror (DM2) is monitored 
by an InGaAs infrared photodetector Djf>. It is used to stabilize the OPO cavity by a servo-loop 
made by a high-voltage amplifler connected to the PZT, controlling the output power of the IR beam 
and stabilizing it during a time ranging from seconds to minutes of continuous operation. The green 
light reflected by DM1 and DM2 is detected by a visible photodetector Dq- 

After the dichroic mirror (DM1) at the OPO output, a flipping mirror (M) allows either the imaging 
of signal and idler far and near flelds on a screen, or the recording of their intensities by two InGaAs 
infrared four-quadrants photodetectors Dl and D2 (ETX505Q from Epitaxx). In this experiment, 
only one quadrant of each photodetector is used. Their quantum efficiency are very close (less than 
1% difference) and equal to 90% ± 5%. In both configurations, polarizing beamsplitters (PBS) 
separate the signal and idler beams. Depending of the orientation of the half wave plate (HWP) 
placed before PBSl, signal and idler can be sent either to Dl or D2. During the measurement process, 
two irises of variable diameter are used to select in the far field a circular region of the output beam. 
Iris II is used to select a narrow circular region on both signal and idler, while iris 12 acts only on a 
single beam, either signal or idler. The lenses LI and L2 are adjusted for each experiment in order 
to project the far field of the beam into the iris plane. The transmitted intensities are recorded and 
their fluctuations are monitored by a dedicated data acquisition system that is described below. 
The data acquisition system used for the noise measurement differs from the usual method of direct 
measurement and subtraction of noise fluctuations (e.g. Ref.[^). The high frequency part of the 
photocurrent of each photodetector is amplifled by a transimpedance amplifler and a broadband 36 
dB amplifler. This signal is then mixed in an electronic demodulator at a frequency /o equal to 
3.5 MHz which lies inside the cavity bandwidth for the infrared modes and above the excess noise 
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frequencies of the output beams. The output of the mixer has an active low-pass filter, working 
at 100 kHz. This signal is registered by a fast Analog-to-Digital acquisition card for measuring the 
noise correlation of the photocurrents. We used two acquisition cards (PCI6110E from National 
Instruments) with four simultaneous measurement channels each and 12 bits for signal conversion. 
The signal is measured with a repetition rate of 200 kHz and registered in the computer. Remaining 
channels of the data acquisition system are used for the measurement of the average value of the 
photocurrent of each detector, as well as the photocurrents of detectors Dq and Dir. Prom the 
stored information of the average value and noise fluctuations of the photocurrent we can calculate 
the noise correlation of the intensity, and compare it to a previous calibration of the shot noise level 
made with a single output of the OPO or with the injected IR light from the Nd:YAG laser. This 
technique allows us to acquire in a very short time interval both the temporal fluctuations of a given 
Pourier component of the photocurrents and of their mean values, and then to postprocess the stored 
data. We can thus determine the different normalized quantities that we have defined earlier in that 
paper. 

C. Experimental procedure 

When the OPO oscillation is stabilized for a given cavity length, we perform two kinds of experi- 
ments. In the first experiment, we record the simultaneous noise distribution of the signal and idler 
beams by using the iris II. We continuously close the iris II in 2s, while acquiring a long series of 
values of Sig, Sii, is, ii, iiR and io (6 series of 400,000 simultaneous values). We analyze the data 
by groups of 10,000 values, calculating for each group the average photocurrent value (i^) (m=l,2) 
and the normalized variances ricorr and na- The corresponding vacuum fluctuation A^E^ was pre- 
viously calibrated, and the electronic noise level is subtracted in the calculation of the variances. 
The transmittance of the iris for each series of 10000 points is deflned as T = where r is the 
ratio r = ^^(i^^f ; fopen corresponds to the initial value of r when the diaphragm is open. The 
normalized intensity difference noise variances ricorr a-nd Ud are then plotted as a function of the 
transmission of the iris. 

The second kind of experiments consists in setting iris I2 on either the signal or the idler beam path 

and taking the unaffected beam as a reference. The procedure is identical to the previous one, except 
that r, for the iris on the signal beam, is deflned as r = -^f^- Closing the iris I2 attenuates only one 
of the beams, and produces a strong imbalance between the photocurrents of the two photodetectors. 
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with a ratio {is/ii) well out of the range 90% to 110%, which was typical in the previous conditions. 
With such a large imbalance, the value of deduced from the experiment does not give any quan- 
titative information on the quantum correlations between the two beams. On the other hand, the 
calculation of the corrected value of the normalized noise ricorr allows one to recover the information 
on the quantum correlation only in the case where the two beams are single mode. Therefore, in this 
second kind of experiments, we only calculate 

COW 3j function of the transmittance T. 



V. SPATIAL DISTRIBUTION OF THE INTENSITY CORRELATION 
A. Experiments with the iris on both signal and idler beams 

As described in the previous subsection, only iris Ji is present. We have realized in this configura- 
tion two series of measurements: in the first one, the OPO oscillation is stabilized for a cavity length 
outside the confocality range, ensuring a TEMqq output mode for signal and idler beams (ref.|T9|, |23[] ) 



; in the second one, the cavity length is set to a shorter value, inside the confocality range, where 
complicated patterns can be observed in the signal and idler beams. While closing continuously the 
iris Ji, we record the transmitted intensity of the signal and idler modes and the normalized intensity 
noise difference variance n^. The different normalized variances are then plotted as a function of 
the transmittance T. In the figures the straight line represents the value of that could be 
calculated with a single mode beam having the same intensity correlation as the whole beam. One 
observes that outside the confocality range, where the signal and idler fields are close to TEMqq, 
the calculated straight line fits well the experimental points. Moreover, for an open diaphragm, 
the normalized intensity difference noise variance is equal to 0.8, demonstrating significant quantum 
intensity correlations between signal and idler. 

In the confocality range, where the signal and idler fields have more complicated, and different, 
transverse variations, one observes that the variation of na is no longer linear with the transmit- 
tance, demonstrating that the emission of the OPO is multimode. Moreover, when the iris is open 
in that confocal configuration, the noise correlation is still below shot noise and close to the value 
obtained for the single mode emission: the quantum intensity correlations for the whole beams are 
preserved even for multimode beams. The increase of noise when the iris is closed shows that the in- 
tensity correlation is stronger in the outer parts of the beams than in the center. Unfortunately, with 
such an experiment it is difficult to describe more precisely the transverse distribution of correlations 
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FIG. 3: spatial distribution of the normalized intensity difference noise Ud as a function of the transmittance 
T with iris Ii on signal and idler beams. Points: values of for a cavity length equal to Lconf = L + 0, 5mm 
(beyond the confocality range)and Lconf = L — 0,4:mm (inside the confocality range). Straight line: Values 
of rid that would be obtained with a single mode beam having the same squeezing when the iris is fully 
open. 

between the signal and idler beams. 

B. Experiments with the iris either on signal or idler beam 

1. Single mode beams 

For a length of cavity L = L^onf + 5, 6mm, where the OPO is supposed far from transverse 
degeneracy, we compared the results of the two kinds of experiments for the values of 
function of the transmittance. The results are plotted in Fig.^. For all those experiments, as the 
beam intensity on the photodetector is reduced, the fluctuations of the photocurrent approaches the 
electronic level. The dispersion in the calculated data increases when we close the iris. Typically, the 
variance of the intensity noise reached the dark noise level for an incident power of 0.2 mW, obtained 
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FIG. 4: Experimental values of ricorr for quasi-TE^Moo signal and idler beams, a) Iris on both beams b)Iris 
on a single beam. Idler (•) and signal(o). Dashed lines: average corrected noise. 

with a typical iris transmittance of 10%. 

One observes that all the experimental points are aligned on a horizontal straight line, with a mean 
value around 70%, which shows that the signal and idler beams are single mode, quantum-correlated, 
beams. The difference in the obtained level of noise for in fig^D (69% and 75%) is certainly due to the 
fact that both experiments are performed at different moments, and therefore on possibly different 
longitudinal modes. The values of ricorr remain stable during a single series of measurements when 
we close the iris. 



2. Beams emitted by an OPO close to confocality 

The length for the cavity is chosen to be very close to the exact confocal length Lconf- In this 



configuration complex structures appear in far and near fields of the signal and idler beams |19]. We 
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performed the same experiments and analysis as in the previous subsection, adding the image of the 
near field and the far field of the beam obtained with the CCD cameras. FigJ^a, b and c display 
respectively the results obtained for L — L^onf = 0.38, —0.37 and —0.62mm. In the configuration of 
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FIG. 5: Experimental values of ricorr for signal and idler beams displaying complex transverse distributions. 
The distance from confocality (AL) is equal respectively to a) 0.38, b) -0.37 , and c) -0.62 mm. Left, iris 
on both beams: (■). Center, iris on a single beam. Idler (•) and signal(o). Right: Far field (FF) and near 
field (NF) images of signal and idler beams. 



Fig.|]a, one observes a nearly gaussian distribution of the idler beam and a non-gaussian distribution 
of the signal beam. When the iris I2 is placed on the idler beam, ricorr does not deviate in a clear way 
from a flat line, except for lower values of the iris transmittance. In this situation, the dispersion in 
the results increases, following the attenuation of the idler field intensity. When we set the iris on the 
signal path, the corrected noise deviates slightly from a constant value for values of the transmittance 
close to 50%. The same behavior can be observed when we perform the experiment with iris 1 on 
both signal and idler. 

In the configuration of Fig.^D, more complex structures appear : a ring pattern can be seen in the 
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outer part of the idler far field image. In the near field, complex structures appear for both signal and 
idler beams. In this situation, when iris Ji is put on both beams, the corrected correlation noise ricorr 
goes up to the shot noise level and stabilizes around this value. It is interesting to observe that the 
level of squeezing obtained is reduced, and the normalized noise for open iris starts at 85%. When 
the iris I2 is put on the signal path and closed, ricorr remains constant: from these observations, one 
can infer that the photons in the signal beams which are correlated to the idler photons are uniformly 
distributed inside the signal beam. On the contrary, when the diaphragm I2 is set on the idler beam, 
the plot of ricorr deviates from a flat line. In that case the results are very close to those obtained 
when the iris Ji is set on both signal and idler, showing that most of the correlated photons in the 
idler are concentrated in the external part of the beam. The center part of the idler is constituted 
essentially of non-correlated photons. 

In the configuration of Fig.^, both signal and idler beams present a faint external ring out of 
the central maximum. In that case, the evolution of the noise when closing iris 1 or closing iris 
2 for signal or idler are very close. Both signal and idler corrected noises ricorr are reduced when 
one closes iris 12, showing that unlike the centers of signal and idler beams, their outer parts are 
quantum correlated. We see that these experiments give interesting indications on the details of 
transverse distribution of the correlations inside the beams. Unfortunately, the great quantity of 
modes oscillating simultaneously in a confocal cavity, which makes a theoretical approach of the 
system very difficult, and the lack of long term stability of the setup, which induces a large dispersion 
on the experimental results, prevent us from a more quantitative comparison between theory and 
experiments, which has been tackled in the case of the multimode VCSEL for example |§. 

VI. CONCLUSION 

The results presented in this paper show that the intensity correlations in a confocal OPO are at 
the same time below the standard quantum limit and not uniformly distributed inside the beams, 
which is a clear evidence that the quantum state generated by such an OPO is a multimode non- 
classical state of light. This experiment showed also that even in the multimode case, the intensity 
fluctuations of the signal and idler beams remain quantum correlated. It would be very interesting to 
find a theoretical explanation, even qualitative, to the fact that, in our experimental configuration, 
the central part of the beams seem to be less quantum correlated than the outer parts. We plan to 
extend this study in the case of an OPO working in a lower degeneracy cavity, where only some few 
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modes are allowed to oscillate simultaneously. The small number of oscillating modes can lead to 
a detailed theoretical description of the system. Another interesting regime is the operation below 
threshold, where the comparison with theory is somehow simpler, but which requires a spatially 
resolved homodyne detection. For a more general point, this demonstration of an intrinsic multi- 
mode emission of signal and idler beams in the quantum regime of the above threshold OPO, opens 
interesting prospects for the study and use of transverse quantum effects in bright beams. 
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